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Overview 


We  propose  a  new  approach  for  parametric  modeling  and 
estimation  of  STAP  data,  based  on  the  2-D  Wold  like 
decomposition  of  random  fields. 

•  The  proposed  parametric  estimation  algorithms  of 
the  interference  components  simplify  and  improve 
existing  STAP  methods. 

•  The  resulting  modeling  and  processing  methods  pro¬ 
vide  new  parametric  tools  to  estimate  and  mitigate 
the  Doppler  ambiguous  clutter. 

•  The  estimation  algorithms  we  propose  enable  the 
estimation  of  the  interference  signals  using  the  ob¬ 
servations  in  only  a  single  range  gate. 

•  The  proposed  method  is  particularly  suitable  for 
non- stationary  clutter  and  jamming  environments. 

•  The  approach  also  provides  a  new  analytical  insight 
into  the  STAP  problem. 


The  STAP  Problem 


The  goal  of  space-time  adaptive  processing  is  to  ma¬ 
nipulate  the  available  data  to  achieve  high  gain  at 
the  target  angle  and  Doppler  and  maximal  mitiga¬ 
tion  along  both  the  jamming  and  clutter  lines. 


Because  the  interference  covariance  matrix  is  un¬ 
known  a  priori,  it  is  typically  estimated  using  sam¬ 
ple  covariances  obtained  from  averaging  over  a  few 
range  gates. 


The  Approach 


We  adopt  the  2-D  Wold  like  decomposition  of  ran¬ 
dom  fields  as  the  parametric  model  of  the  observed 
data. 

Employing  this  model,  we  derive  computationally 
efficient  algorithms  useful  for  parametrically  esti¬ 
mating  both  the  jamming  and  clutter  fields. 

Having  estimated  the  interference  terms  parametric 
models,  their  covariance  matrix  can  be  evaluated 
based  on  the  estimated  parameters. 

Once  the  parametric  models  of  the  interference  com¬ 
ponents  have  been  estimated,  several  alternatives 
for  mitigating  the  interference  are  available. 

•  Parametric  fully  adaptive  processing. 

•  Parametric  partially  adaptive  processing  -  only 
the  spectral  support  parameters  of  the  interfer¬ 
ence  components  need  to  be  estimated. 


The  2-D  Wold-Like  Decomposition 

Theorem:  Let  {y(n,  m),  (n,  m)  G  .2^^}  be  a  regular  ran¬ 
dom  field.  Then  {^(n,  m)}  can  be  represented  uniquely 
by  the  orthogonal  decomposition: 

y(n,  m)  =  rc(n,  m)  +  u(n,  m)  , 


w{n^m)  =  a{k^l)u{n  —  k^m  —  1)  ^ 


u{n^m)  is  the  white  innovations  field  of  y(n,  m). 

The  field  {rc(n,  m)}  is  purely-indeterministic  and 
regular. 


The  field  {u(n,  m)}  is  deterministic. 


The  Prediction  Support 


RNSHP  support;  example  with  o  =  2  and  [3=1. 


The  Deterministic  Component  Decomposition 


v{n,  m)  =  h{n,  m)  +  ^  (n,  m) 

GO 

where  m)}  is  a  harmonic  random  field  and 

{e(Q;,/3)(n,  m)}  is  an  evanescent  field  . 

O  is  a  set  of  pairs  (o,  /?)  of  coprime  integers. 


All  decomposition  components  are  mutually  orthogonal. 


Modeling  the  Evanescent  Field 


For  any  (o,  /?)  G  O 
e{a,0){n,m)  = 


i 

The  1-D  purely-indeterministic,  complex  valued  pro¬ 
cesses  {s\^'^\na  —  m/3)}  and  —  m/3)},  are 

zero- mean  and  mutually  orthogonal  for  all  i  ^  j. 


The  “spectral  density  function”  of  each  evanescent  field 
has  the  form  of  a  sum  of  1-D  delta  functions  which  are 
supported  on  lines  of  rational  slope  in  the  2-D  spectral 
domain. 


The  Spectrum  of  Evanescent  Fields  —  Example 


DFT  of  an  observed  field  with  two  evanescent  components. 


The  Harmonic  Component 


h{n,  m)  =  Y^ 


P 


The  Cp  ’s  are  mutually  orthogonal  random  variables, 
E\Cp\‘^  =  cTp,  and  {ujp^Up)  are  the  spatial  frequencies  of 
the  p  —  th  harmonic. 


STAP  and  the  2-D  Wold  Decomposition 


The  target  signal  is  modeled  as  a  random  amplitude 
complex  exponential  -  in  the  space-time  domain  the 
target  model  is  that  of  a  2-D  harmonic  component. 


The  purely-indeterministic  component  of  the  space- 
time  field  is  the  sum  of  a  white  noise  field  due  to  the 
internally  generated  receiver  amplifier  noise,  and  a 
colored  noise  field  due  to  the  sky  noise  contribution. 

The  interference  components  (jammers  and  clutter) 
are  evanescent  random  fields: 

•  In  the  space-time  domain  each  jammer  is  mod¬ 
eled  as  an  evanescent  component  with  {a,  (3)  = 
(1,0)  whose  1-D  modulating  process  is  white 
noise. 

•  Clutter  from  all  angles  lies  in  a  “clutter  ridge” , 
supported  on  a  diagonal  line  (that  generally 
wraps  around  in  Doppler)  in  the  angle-Doppler 
domain,  with  (cr,/?)  such  that  tan/3/(a  corre¬ 
sponds  to  the  slope  of  the  clutter  ridge. 


Setting  the  Estimation  Problem 


Let  {^(n,  m)},  (n,  m)  G  D  where  D  =  {(i,j)|0  <  i  < 
S  —  IjO  <  j  <T  —  1}  be  the  observed  random  field. 


Define 


y=  b(0,0),...,y(0,T-  l),y{l,0),...,y{l,T  -  1) 


P  • 


{a,f3) 


P  • 


(o:,/3} 


Setting  the  Estimation  Problem 


Let 


^(a,/3)  _  1-^(0:,^) 


s;  -  ^\sr^>(0),s\ 


{a, (3) 


{a, (3) 


((5  —  1)0)  .  .  . 


(a,/3) 


-/3)...,4“’^)(a-(r-l)/?)... 


be  the  vector  of  the  observed  samples  from  the  1-D  mod¬ 
ulating  process  of  the  evanescent  field 


Define 


y{oL,fi)  =  [  0,  a,  . . . ,  (T  —  l)a, 

/3,  j3  -\~  o,  •  •  •  5  P 

(5-1)/?,  (5- !)/?  +  «,...  (5-l)/?+(r-l)a]^ 


Then 


(oi,(3) 


q;^+/3^ 


r’' 


Properties  of  the  Evanescent  Field 


Whenever  na  —  mp  =  ka  —  ip  for  some  integers  n,  m,  fc,  I 
such  that  {)  <  n^k  <  S  —  1  and  0  <  m, ^  <  T  —  1,  the 

same  element  of  appears  more  than  once  in  the 

vector. 


Lemma:  For  a  rectangular  observed  field  of  dimensions 
S  X  T  the  number  of  distinct  samples  from  the  random 

process  that  are  found  in  the  observed  field  is 

{S  —  l)|o|  +  (T  —  l)\P\  +  1  —  (|o|  —  1)(|/3|  —  1) 


The  concentrated  version^  of  is  a  (5  —  1)101  + 

(T  —  1)|/3|  +  1  —  (|o|  —  1)(|/3|  —  1)  column  vector  of  non¬ 
repeating  samples  of  the  process 


poc,(3) 


where  is  rectangular  matrix  of  zeros  and  ones 

which  replicates  rows  of 


The  Covariance  Matrix  of  the  Field 


The  covariance  matrix  which  characterizes 

L 

{^(«,/3)}  .g  terms  of  the  concentrated  version 

vector 


The  covariance  matrix  of  the  evanescent  component  is 
given  by 


-p(a,/3) 

^  i 


{oc,(3)  f  i{oc,(3)xh- 


Since  the  evanescent  components  are  mutually 

orthogonal,  and  since  all  the  evanescent  components  are 
orthogonal  to  the  purely-indeterministic  component 


r  =  Tp/  +  5]  5] 

(a,l3)eO  i=l 


Estimating  the  Interference  Parametric  Model 


•  Estimate  the  (o,  /?)  pair. 

•  For  a  fixed  c  =  na  —  m/3  (i.e.,  along  a  line  on  the 
sampling  grid),  the  samples  of  the  evanescent  com¬ 
ponent  are  the  samples  of  1-D  constant  amplitude 

harmonic  signal,  whose  frequency  is 


E  y{n,m)exp{-j2Tr 


na.—ml3=c 


OL^  -h 


(n/3  -h  ma 


Ns  denotes  the  number  of  the  observed  field  samples  that 
satisfy  the  relation  na  —  mf3  =  c. 

•  Estimate  the  parameters  of  the  1-D  AR  (or  any 
other)  model  of 


Fully  Adaptive  Parametric  Processing 


Having  estimated  the  parametric  models  of  the  noise  and 
interference  components  of  the  field,  the  estimated  pa¬ 
rameters  are  substituted  into  the  parametric  expression 
of  the  covariance  matrix  to  obtain  an  estimate  of  the 
interference-plus- noise  covariance  matrix  F. 

Weights  in  the  optimal  space-time  filter  are  given  by 

w  =  T~^vt 


The  test  statistic  is 


z{vo^  id)  =  {vo^  id)y 
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Partially  Adaptive  Parametric  Processing 


After  demodulation  using  the  estimated  o,  /?, 

p(Q;,/3)  _ 

%  %  Z  V  2'  / 


Construct  the  following  orthogonal  projection  matrix 


(«,/3)  ^  ,(« 


m 


It  is  easily  verified  (by  substitution)  that  is  an 

orthogonal  projection  onto  the  range  space  of  since 

^  i  ^  ^  i  -‘-i 


The  projection  matrix  onto  the  subspace  orthogonal  to 
the  clutter  space  is  given  by 

•  is  a  sparse  matrix  of  zeros  and  ones,  whose 
structure  is  a  function  of  a  and  p  only.  Hence, 

jg  easily  computed. 

•  Only  a,  /?,  need  to  be  estimated  ! 


